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Technical  Progress  Report 


1.  Introduction 

The  goal  of  this  research  is  to  increase  our  understanding  of  the  impact  of  the  ocean 
and  seafloor  environmental  variability  on  deep-water  (long-range)  ocean  acoustic 
propagation  and  to  develop  methodologies  for  including  this  in  acoustic  models. 
Experimental  analysis  is  combined  with  model  development  to  isolate  specific 
physics  and  improve  our  understanding.  During  the  past  few  years,  the  physics 
effects  studied  have  been  three-dimensional  propagation  on  global  scales,  deep 
water  ambient  noise,  under-ice  scattering,  bathymetric  diffraction  and  the 
application  of  the  ocean  acoustic  Parabolic  Equation  to  infrasound. 


2.  Tasks 


a.  Task  1:  Basin  Scale  Acoustics  and  CTBTO  Data  Analysis 

Scenario 

The  scenario  chosen  to  evaluate  the  impact  of  mesoscale  variability  on  the  arrival 
angle  of  long-range  signals  is  a  seismic  event  on  the  Kerguelen  Plateau  (-53°S  71°E) 
in  the  southern  ocean.  This  region  of  the  world,  which  includes  Heard  Island,  Crozet 
Island  and  the  Kerguelen  Plateau  has  historical  significance  in  the  long-range 
underwater  community  (Perth-Bermuda,  HIFT)  and  is  the  position  of  the  IMS 
hydro-acoustic  station  HA04.  As  demonstrated  in  the  Heard  Island  Feasibility  Test, 
sound  sources  here  can  be  detected  from  the  Southern  Atlantic  to  the  Pacific.  The 
receiver  position  is  taken  be  approximately  at  the  IMS  HA10S  location  (Ascension 
Island  southern  station).  The  hydroacoustic  station  at  Ascension,  as  can  be  seen  in 
Figure  1,  has  acoustic  visibility  to  the  South  Atlantic,  South  Indian  and  South  Pacific 
and  has  a  large  number  of  hydroacoustic  signals  generated  by  seismic  events. 

The  Kerguelen-Ascension  path  was  also  chosen  because  of  the  relative  minimal 
bathymetric  variability  along  the  geodesic.  For  low-frequency  acoustic  signals, 
bathymetric  variation  leads  to  refraction1  and  diffraction2  and  this  complicates  the 
evaluation  of  the  impact  of  mesoscale  variability.  Mesoscale  variability  refracts 
sound  from  horizontal  gradients,  primarily  ocean  fronts  and  eddies.  This  scenario 
has  acoustic  energy  crossing  the  Antarctic  Circumpolar  Convergence,  the  strong 
boundary  between  the  Southern  Ocean  and  the  Indian  Ocean.  This  path  also 
transects  the  Agulhas  Retroflection,  a  dynamic  region  of  eddies  that  spin  off 
interaction  of  the  Agulhas  current  with  the  coast  of  South  Africa3.  The  deflection  of 
acoustic  energy  due  to  interaction  with  these  eddies  was  evaluated  by  Munk  to  see  if 
they  could  explain  the  Perth-Bermuda  results4.  He  found  them  to  be  too  weak  for 
the  geometry  of  Perth-Bermuda. 
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The  time  and  range-dependent  ocean  used  for  this  simulation  study  is  the  eddy¬ 
resolving  ECC02  (Estimating  the  Circulation  and  Climate  of  the  Ocean,  Phase  II 
(ECC02):  High-Resolution  Global-Ocean  and  Sea-Ice  Data  Synthesis)  model  re¬ 
analysis  for  the  years  1992  and  1993.  The  ECC02  model  is  a  state  estimation  based 
upon  data  syntheses  obtained  by  least  squares  fitting  of  the  global  ocean  and  sea-ice 
configuration  of  the  Massachusetts  Institute  of  Technology  general  circulation 
model  (MITgcm)  to  the  available  satellite  and  in-situ  data5.  The  ECC02  product  is 
the  temperature  and  salinity  vs.  depth  (50  depth  layers)  on  a  %  degree  grid, 
produced  every  3  days.  The  sea-surface  temperature  (SST)  for  ECC02  model  for 
January  1, 1992  is  shown  in  Figure  2.  Regions  of  strong  mesoscale  activity  (eddies) 
are  visible  of  the  Northeast  coast  of  the  United  States  (the  Gulf  Stream)  and  Japan 
(the  Kuroshio)  and  south  west  of  South  Africa  (the  Agulhas  Retroflection). 

Ecco2  Model  (19920105) 


Figure  2.  ECC02  sea-surface  temperature  for  19920101. 

Modeling  approach 

The  Peregrine  model  is  used  to  compute  the  acoustic  field,  including  out-of-plane 
propagation  with  the  ECC02  reanalysis  model  every  3-days  from  1992  to  1993.  In 
order  to  reduce  the  computational  cost  for  long-range  propagation  in  Peregrine,  the 
computational  grid  can  be  confined  to  a  swath  around  a  geodesic  bearing.  The 
center  bearing  is  268°  and  the  swath  width  is  1200  km.  The  source  depth  is  730m 
and  the  maximum  range  computed  is  9100  km.  The  bathymetry  is  taken  from  the 
ETOPOl  dataset.  For  the  2Hz  computation,  where  even  the  deep-water  bathymetry 
did  show  an  impact  on  the  field,  a  flat  seafloor  of  3000m  depth  was  imposed  to 
remove  this  effect.  The  seafloor  is  set  to  be  an  acoustically  deep  layer  of  medium  silt 
(sediment  grain-size  6),  and  is  not  expected  to  play  a  role  in  this  study. 
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Kerguellen  Plateau  to  Ascension  3D  Propagation 


Figure  3.  8Hz  depth  averaged  Transmission  Loss  (TL)  for  Kerguelen  source  showing 
propagation  tube. 

The  depth  averaged  transmission  loss  in  Figure  3  shows  the  diffractive  filling  of  the 
shadows  behind  Crozet  Island  and  Agulhas  shelf,  as  well  as  the  open  water  path 
from  Kerguelen  to  Ascension  (the  center  of  the  final  range  swath).  The  field  was 
computed  with  4  points  per  wavelength  in  the  horizontal  and  4  Pade  terms  in  the 
Split-Step  Pade6  expansion  in  the  horizontal.  The  field  was  computed  every  3  days 
for  2,  4  and  8  Hz  and  every  month  for  16  Hz. 

In  order  to  estimate  the  back-azimuth,  a  2D  planar  array  centered  at  8°  53’  S  and  14° 
37’  30"  W  was  used.  The  complex  pressure  field  was  computed  on  a  360  x  312  grid 
with  spacings  of  30.9m  and  35.6m  for  longitude  and  latitude,  respectively.  The 
receiver  depth  was  834m,  the  reported  depth  of  HA10S.  Narrow  band  plane  wave 
beamforming,  with  Hann  window  shading,  was  performed  on  each  arrival.  The 
geodesic  back-azimuth,  using  the  WGS’84  globe,  between  the  source  and  array 
center  is  142.498°.  For  all  frequencies  computed,  the  array  spacing  is  greater  than  a 
quarter-wavelength,  but  grating  lobes  are  handled  by  only  beamforming  around  the 
expected  azimuth  of  the  acoustic  energy. 
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Results 

Prior  to  evaluating  the  mesoscale  refraction  of  the  eddy-resolving  ECC02  model,  the 
refraction  induced  by  large-scale  oceanographic  features,  such  as  crossing  the 
Antarctic  Circumpolar  Front  is  examined.  The  World  Ocean  Atlas  20097 
temperature  and  salinity  fields  for  the  four  seasons  was  used.  The  results,  for  2,  4,  8 
and  16  Hz,  are  plotted  in  Figure  4,  along  with  the  geodesic  (in  cyan).  There  is  very 
little  seasonal  dependence  observed  and  all  of  the  back-azimuths  are  within  0.25°  of 
the  geodesic.  This  minimal  deviation  from  the  geodesic  may  result  from  the 
inherent  smoothing  of  gradients  when  computing  a  climatology  database. 


World  Ocean  Atlas  Backazimuth  2(k),  4(b),  8(g),  16(r)  Hz  /  Geodesic  (c) 


Figure  4.  Centroid  ofbeamformed  arrival  vs  season  using  the  World  Ocean  Atlas  for 
each  season  (Winter,  Spring,  Summer  and  Fall).  The  frequencies  computed  are  2Hz 
(black),  4Hz  (blue),  8Hz  (green)  and  16  Hz  (red).  The  un-refracted geodesic  is  the 
thick  line  in  cyan  and  142.498°. 

The  acoustic  field  for  2,  4  and  8  Hz  was  run  through  the  ECC02  model  from  January 
3, 1992  through  December  28, 1993,  sampled  in  three  day  increments,  yielding  a 
total  of  243  sample  time  series.  Due  to  the  frequency-cubed  cost  of  the  3D  PE 
computation,  the  16  Hz  runs  were  computed  monthly  for  the  same  time  period. 

The  beam  time-series,  at  the  output  of  the  beamformer  for  the  4Hz  signal  over  the 
back-azimuths  from  135  to  150°  is  shown  in  Figure  4.  The  received  beam  energy  is 
generally  unimodal,  with  a  smooth  distribution  around  a  center  beam,  varies  in 
arrival  angle  and  is  dynamic  in  level.  The  dynamic  range  of  the  plot  is  30  dB. 
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Ecco2  Kerguellen-Ascension  3D  4Hz 


2  4  6  8  10  12  14  16  18  20  22 

Time  From  19920105  (Months) 


Figure  5.  Beam  response  for  4Hz  propagation  over  the  1992-1993  Ecco2  model. 

The  beamformer  output  was  collapsed  to  a  single  back-azimuth  value  by  computing 
the  centroid  of  the  beam  power  over  the  searched  arrival  angle.  The  time-series  of 
the  back-azimuth  for  the  set  of  model  simulations  is  shown  Figure  6,  along  with  the 
geodesic  back-azimuth.  This  simulation  is  not  finished  to  the  point  of  being 
publishable.  There  is  substantially  more  refraction  observed  in  the  ECC02  result 
than  in  the  WOA09  result,  indicating  the  combined  impact  of  the  mesoscale  eddies 
and  sharper  front  definition  is  important  to  the  acoustic  propagation  path.  There  is 
a  long  time-scale  oscillation  on  the  order  of  a  120-day  period,  ostensibly  related  to 
seasons,  but  it  does  not  repeat  for  each  season  of  1992  to  the  corresponding  season 
of  1993.  This  oscillation  leads  to  the  range  in  the  back-azimuths  going  from  a 
minimum  near  141.5°  to  a  maximum  of  143.1°.  The  back-azimuth  as  a  function  of 
frequency  does  show  coherent  behavior,  at  least  within  the  observed  0.3°  small  time 
scale  variability.  There  is  a  consistent  3-6  day  oscillation  with  a  peak  magnitude  of 
0.3°  in  the  4  and  8Hz  data  which  could  be  random  noise,  or  could  be  due  to  small 
scale  motions  of  the  eddies  and  the  bathymetry.  Note  the  2Hz  signal  does  not  have 
this  oscillation.  Recall  that  it  was  computed  using  a  flat  seafloor  to  remove 
dominant  bathymetry  effects. 
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This  needs  to  be  finished 


Kerguellen-Ascension  3D  2(k),4(b),8(g).16(r)  Hz 


Figure  6.  Centroid  of  arrival  energy  for  2-16  Hz  over  1992-1993  Ecco2  model. 

b.  Task  2:  NPAL  PhilSealO  Data  Analysis  and  Matched  Field 
Processing 

PhilSea09  analysis  involved  noise  modeling  and  analysis  of  the  vertical  directivity 
from  the  PS09  Vertical  Line  Array.  PS09  was  chosen  because  this  VLA  was  spaced 
at  F/2  at  250Hz  and  is  therefore  beamforming  capable,  covering  the  conjugate 
depth. 

An  ambient  noise  model  was  built  for  predicting  the  vertical  noise  component  of  the 
ambient  noise.  The  model  used  Peregrine  for  hydrophone  to  A,  /4  depth  below  the 
surface  (for  a  noise  sheet  and  for  ships)  to  the  edge  of  the  Philippine  Sea  basin. 

Wind  speed  was  computed  by  adding  each  surface  patch  to  the  array  cross-spectral 
density  (corresponding  to  the  Kuperman-Ingenito  surface  sheet  of  independent 
sources).  Shipping  was  input  via  a  realization  of  the  HITS  model  for  the  Philippine 
Sea.  The  results  are  dependent  upon  the  seafloor.  For  a  soft  sediment,  the  modeled 
vertical  noise  directivity  at  the  axis  and  at  the  conjugate  depth  are  shown  in  Figure 
7. 
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Beam  Noise  on  Axial  and  Deep  Subarrays: 

Soft  Bottom 

5-50Hz,  Wind  +  Distant  HVMS,  Axial  Array  BF,  Soft  Bottom  5-200Hz,  Wind  +  Distant  HVMS,  Deep  Array  BF,  Soft  Bottom 


•  No  noise  notch  below  conjugate  region 

•  Moderate  noise  notch  at  higher  freqs  at  axis 


Figure  7.  Vertical  Noise  Directivity  as  a  function  of  frequency  for  an  axial  (left)  and 
deep  (right)  vertical  line  array  for  frequencies  from  5-50Hz  ( axial  array)  and  5-200Hz. 
-  Which  are  the  unaliased  beamforming  capable  frequencies  of  the  deployed  arrays. 

We  now  compare  the  deep  VLA  noise  vertical  directivity  of  the  measurement  from 
those  modeled  using  a  hard  seafloor  and  a  soft  seafloor.  The  result  is  shown  in 
Figure  8,  indicating  the  seafloor  is  best  modeled  as  hard. 


OASIS,  INC. 


9 


Report  No.  QSR-14C0 172-Ocean  Acoustics-063015 


Model/Data  Comparison:  Beams 


1 50Hz,  Wind  +  Dist,  Hard  Bottom  1 50Hz,  Wind  +  Dist,  Soft  Bottom 


Figure  8.  Vertical  Noise  Directivity  for  150Hz  using  a  hard  seafloor  (left),  soft  (middle) 
and  from  the  observations.  The  lack  of  a  clear  critical  angle  indicates  a  harder 
sediment. 

The  deep  VLA  data  was  processed  for  24  hours.  The  band  averaged,  Hann- 
windowed  BTR  (Bearing  Time  Record)  of  vertical  angle  vs.  time  is  shown  in  Figure 
8. 
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angle  (degrees  from  horizontal) 


JD1  13  NPAL  DVLA  25-75  Hz  Noise  Directionality 


1  2  3  4  5  6  7  8  9  1011 

Hours  from  031650Z 

JD1  13  NPAL  DVLA  75-125  Hz  Noise  Directionality 


1  234  56789  lO  11 

Hours  from  031650Z 


Figure  9.  Bearing  Time  Records  for  25-75  Hz  and  75-125Hzfor  the  DVLA  on  Julian  Day 
113.  TheCPAoftheR/VRevelleisseen  at  hour  2:30. 

The  CPA  of  the  surface  ship  is  evident  at  hour  2,  with  the  associated  high  angle 
energy,  particularly  at  low  frequency.  Note  the  energy,  particularly  in  the  upper 
band  during  the  quiet  periods  that  is  coming  in  from  the  above  the  array  (positive 
vertical  angles).  This  is  wind  energy  from  above.  It  is  not  evidently  coming  from 
below,  a  physical  phenomenon  Harrison  attributed  to  the  bottom  energy  having  1 
more  bottom  bounce.  This  leads  to  the  technique  of  estimating  the  reflection 
coefficient  by  subtracting  the  downward  energy  from  the  upward  energy  (for  each 
corresponding  beam).  The  result  of  this  operation  is  the  estimated  reflection 
coefficient  vs.  frequency  and  is  shown  in  Figure  10. 
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Ambient  Noise  based  Reflection 
Coefficient  Estimation 


Reflection  Coefficient  Estimate  from  PhilSea09  NPAL  DVLA  Data 


25  50  75  100  125 

Frequency  (Hz) 


Figure  10.  Estimated  reflection  coefficient  as  a  function  of  frequency  by  taking  the 
difference  of  downgoing  and  upgoing  energy. 


3.  Future  Plans 

For  the  next  quarter,  the  focus  of  the  work  will  be  on  examining  measurements  and 
models  of  the  3-dimensional  diffraction  induced  by  bathymetric  scattering. 
Observations  have  been  made  from  the  CTBTO  hydroacoustic  stations  of  seismic 
events  that  are  in  the  acoustic  shadow  due  to  island  (or  continent)  blockage.  The  3D 
Peregrine  model  will  be  applied  to  these  examples  demonstrating  that  diffraction 
can  explain  the  observations  of  hydroacoustic  signals  in  the  deep  shadow. 

Theoretical  work,  following  Munk's  approach,  will  be  done  to  compute  the  expected 
energy  observed  due  to  diffraction  from  an  edge. 

The  impact  of  bathymetric  diffraction  on  the  global  coverage  of  the  CTBTO  network 
will  be  evaluated.  Early  computations  indicate  that  the  filling  in  from  bathymetric 
diffraction  can  be  on  the  order  of  3%  of  the  earth’s  globe,  for  a  single  station,  which 
corresponds  to  15  million  kmA2. 
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4.  Publications  and  Peer  Interactions 

Dr.  Heaney  presented  at  the  ASA  in  Pittsburgh  and  at  the  Conference  on  Ocean  Noise 
in  Sevilla,  Spain.  He  met  with  Dr.  Laslo  Evers  (KNMI)  and  Dr.  Michael  Ainsle,  during 
a  visit  to  the  Netherlands.  Dr.  Heaney  visited  the  CTBTO  (Vienna)  and  met  with  Dr. 
Mario  Zampolli.  He  wrote  a  conference  paper  and  presented  on  global  acoustic 
propagation  (including  on  Europa,  a  small  moon  of  Jupiter)  at  the  International 
Conference  of  Sound  and  Vibration  (ICSV)  in  Florence,  Italy. 

Dr.  Heaney  has  written  a  Mesoscale  Refraction  paper  for  JASA  and  it  will  be 
submitted  in  the  next  few  weeks. 
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Financial  Progress  Report:  Period  ending  June  30,  2015 

OASIS,  INC. 

JOB  STATUS  REPORT  6/30/201 5 


1172  DEEP  WATER  ACOUSTICS  N00014-1 14-C-0172 

POP:  9/27/13-3/6/16 


CONTRACT  VALUE 

Contract  Value 

Cost 

$368,935 

Fee 

$27,048 

Total 

$395,983 

Funding  Value: 

$215,404 

$15,791 

$231,195 

Remaining  to  Fund:  $153,531 

CUM ULATIVE  SPENDING  WITH  COM M ITM  ENTS 

$1 1 ,257 

$164,788 

DIRECT 

OH 

MH 

TOTL  COST 

FEE 

TOTAL 

ACTUAL 

OASIS 

$84,368 

$62,700 

$1,842.00 

$148,910 

$11,168 

$160,078 

COMMITTED 

$0 

$0 

$0 

$0 

$0 

$0 

$84,368 

$62,700 

$1 ,842 

$148,910 

$11,168 

$160,078 

TOTAL  REMAINING  TO  SPEND: 

$71,117 
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